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QQ . ABSTRACT 

We present the results of a spectroscopic monitoring campaign of the OB-star companions to the eclipsing X-ray pulsars SMC X-1, 
LMC X-4 and Cen X-3. High-resolution optical spectra obtained with UVES on the ESO Very Large Telescope are used to determine 
, the radial-velocity orbit of the OB (super)giants with high precision. The excellent quality of the spectra provides the opportunity to 

measure the radial-velocity curve based on individual lines, and to study the effect of possible distortions of the line profiles due to 
e.g. X-ray heating on the derived radial-velocity amplitude. Several spectral lines show intrinsic variations with orbital phase. The 
magnitude of these variations depends on line strength, and thus provides a criterion to select lines that do not suffer from distortions. 
The undistorted lines show a larger radial-velocity amplitude than the distorted lines, consistent with model predictions. Application 
^ ' of our line-selection criteria results in a mean radial- velocity amplitude A'opi of 20.2 ± 1.1, 35.1 ± 1.5, and 27.5 ±2.3 km s"' (Icr errors), 

' for the OB companion to SMC X-1, LMC X-4 and Cen X-3, respectively. Adding information on the projected rotational velocity 

of the OB companion (derived from our spectra), the duration of X-ray eclipse and orbital parameters of the X-ray pulsar (obtained 
from literature), we aiTive at a neutron star mass of 1.06+j; 1^, 1.25+j; |i and 1.34+j; ]^ Mq for SMC X-1, LMC X-4 and Cen X-3, 
^ ' respectively. The mass of SMC X-1 is near the minimum mass (~1 Mq) expected for a neutron star produced in a supernova. We 

' discuss the implications of the measured mass distribution on the neutron-star formation mechanism, in relation to the evolutionary 

\ history of the massive binaries. 
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^ 1 . Introduction mass: 2.1 + 0.2 M© (Nice et aL 2005). Both results are in favor 

i-^ of a stiff EOS (see also Srinivasan 2001). Neutron stars also have 

' o?r°" ^'^'P''^' 'T'"''"^ f ^ T-T^ T a minimum mass limit. The minimum stable neutron-star mass 

> 8 Mo) with a central density that can be as high as 5 to 10 ^^^^^ q ^ ^ ^^^^^^^ ^ ^^^^ ^^^^-^^-^ minimum stems from 

times the density of an atomic nucleus. The global structure of a ^ ^^^^^^^ ^^-^ ^ supernova. Lepton-rich proto neu- 

neutron star depends on the equation of state (EOS) under these ^^^^ ^^^^^ ^^^^^^^^^ ^^^-^ ^^^^^^ j^^^ ^j^^^^ ^^^^^ ^ ^ 

^ . extreme conditions, i.e. the relation between pressure and den- (Lattimer & Prakash 2004; Haensel et al. 2002). 

■ - ' sity in the neutron star interior (e.g. Lattimer & Prakash 2004). . , • • , 

^. r-.^e J- 1 f »u » » J Another issue IS the neutron-star mass distribution: the de- 

Given an EOS, a mass-radius relation for the neutron star and a ., , .... . , 

corresponding maximum neutron-star mass can be derived. The [^'^^"^ supernova mass ejection mechariism accompanying the 

"stiffness" of the EOS depends e.g. on how many bosons are formation of the neutron star is not uiiderstood, but it is likely 

present in matter of such a high density. As bosons do not con- "^^"^ neutrinos that are produced during the forma- 

tribute to the fermi pressure, their presence will tend to "soften" ^1°" f *^ (P™t°-) neutron star in the centre of the collapsing 

the EOS. For a soft EOS, the maximum neutron-star mass will ^" important role (e.g. Burrows 2000). Timmes et al. 

be low (e.g. < 1.55 Mq for the EOS appHed by Brown & Bethe (1996 present mode calculations from which they conclude that 

(1994)); for a higher mass, the object would collapse into a black ^ype II supernovae (massive, single stars) will give a bimodal 

j^pjg neutron-star mass distribution, with peaks at 1.28 and 1.73 Mq, 

r^, , . f ^ . ■ while Type lb supernovae (such as produced by stars in binaries. 

The accurate measurement of neutron-star masses is there- , ■ , ■ , r , ■ , s . 

f ■ ^ J » J- e e » which are Stripped of their envelopes) Will produce neutron Stars 

fore important for our understanding of the EOS of matter at ... , . , r 

, ■ ... ^ »i »u » • » » within a small range around 1.32 M(T). Despite the fact that It IS 

supra-nuclear densities. Currently, the most massive neutron star ... , ° . V ,,, , , , 

■ V u- • V 1 17 1 V wr. • .1 in a binary, the massive neutron star in Vela X-1 may belong to 

in an X-ray binary is the X-ray pulsar Vela X-1 (Barziv et al. , , , ■ , • ,• ., • 

mm <^ 11 1 1 onmN fi o<;,ai<;a^ t-u the second peak in this mass distribution. 

2001; Quaintrell et al. 2003) with amass of 1.86 + 0.16 Mq. The ^ 

milUsecond radio pulsar J075 1 + 1 807 likely has an even higher Neutron stars are detected either as radio pulsars, single or in 

a binary with a white dwarf or neutron star companion, or as X- 

Send offprint requests to: A. van der Meer ray sources in binaries with a (normal) low-mass (LMXB) or a 

* Based on observations obtained at the European Southern high-mass companion star (HMXB). Presendy, all accurate mass 

Observatory at Paranal, Chile (ESO program 68.D-0568) determinations have been for neutron stars that were almost cer- 
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Table 1. Orbital parameters of SMC X-1, LMC X-4 and Cen X-3 obtained from Wojdowski et al. (1998), Levine et al. (2000) and Nagase et al. 
(1992), respectively. To is the mid-eclipse time and corresponds to orbital phase = 0.0. The eccentricity is obtained from Bildsten et al. (1997) 
and references therein. All the reported errors are Icr values unless stated otherwise. For SMC X-1 the range of semi-eclipse angle 9^ includes 
observations obtained by Primini et al. (1976), Bonnet-Bidaud & Van der Klis (1981) and Schreier et al. (1972a); for LMC X-4 those of Li et al. 
(1978), White (1978) and Pietsch et al. (1985); for Cen X-3 we Ust the result of Clark et al. (1988). 





SMC X-1 


LMC X-4 


Cen X-3 


To (MJD) 


42836.18278(20) 


51110.86579(10) 


40958.35(1) 


Porb (days) 


3.89229090(43) 


1.40839776(26) 


2.08713845(5) 


PoJPorb (yr->) 


-3.353(14) X 10-« 


-9.8(7) X 10-^ 


-1.738(4) X 10-« 


^spin (sec) 


0.708 


13.5 


4.82 


fix sin / (It- sec) 


53.4876(4) 


26.343(16) 


39.56(7) 


e 


< 4 X 10-5(2o-) 


0.006(2) 


< 1.6 X 10-3(90%) 


(deg) 


26 - 30.5 


25-29 


32.9 ± 1.4 


OB companion 


BO lb 


08 III 


06.5 II-III 



tainly formed in Type lb supemovae and that have accreted Uttle 
since. Exceptions are J1909-3744, a pulsar (-1- white dwarf) with 
a mass of 1 .438 ±0.024 M© (Jacoby et al. 2005), and the massive 
neutron star in J075 1 -i- 1 807, which may have originated from an 
LMXB. The most accurate masses have been derived for the bi- 
nary radio pulsars. Until recently, all of these were consistent 
with a small mass range near L35 M© (Thorsett & Chakrabarty 
1999). 

We focus here on the initially most massive systems, which 
consist of a massive OB supergiant and a neutron star or a black 
hole (Kaper 2001; Kaper & Van der Meer 2005). The main 
motivation to concentrate on these systems is that they are the 
most likely hosts of massive neutron stars. About a dozen of 
these systems are known'; five of them contain an eclipsing X- 
ray pulsar. The masses of all but one (Vela X-1) are consistent 
(within their errors) with a value of about 1.4 M©. However, 
most spectroscopic observations used for these mass determina- 
tions were carried out more than 20 years ago, before the ad- 
vent of sensitive CCD detectors and 8m-class telescopes, which 
allow high-resolution spectroscopy of the optical companions. 
The uncertainties in the earlier radial velocity measurements (see 
Van Kerkwijk et al. 1995a) are too large to measure a significant 
spread in mass among these neutron stars, if present. 

In this paper we present new, more accurate determinations 
of the mass of the neutron star in three of these systems, i.e. 
SMC X-1, LMC X-4, and Cen X-3 using the high-resolution 
Ultraviolet and Visual Echelle Spectrograph (UVES) on the ESO 
Very Large Telescope (VLT). These systems are in a phase of 
Roche-lobe overflow (Savonije 1978, 1983), have well deter- 
mined, circular orbits (Porb of a few days), and an optical coun- 
terpart of y ^ 14 mag, i.e. well within reach of VLT/UVES. 

In Sect. 2 we introduce the three HMXBs. In Sect. 3 we de- 
scribe the acquired observations and data reduction procedure. In 
Sect. 4 we present the spectral analysis and the resulting radial- 
velocity curves. In Sect. 5 we evaluate the measured radial- 
velocity amplitudes and derive the mass of the neutron star in 
these three systems. In Sect. 6 we summarise our conclusions 
and in Sect. 7 we compare these to the predictions of supernova 
models. 



' Recently, several new hard X-ray sources have been detected with 
INTEGRAL that show the characteristics of a heavily obscured HMXB 
with an OB-supergiant companion, called supergiant fast X-ray tran- 
sients (Negueruela et al. 2006; Lutovinov et al. 2005) 



2. Eclipsing high-mass X-ray binaries; review of 
earlier work 

Five high-mass X-ray binaries are known to host an eclipsing 
X-ray pulsar: Vela X-1, 4U 1538-52, SMC X-1, LMC X-4 
and Cen X-3. The eclipse provides an important constraint on 
the orbital inclination i, an essential parameter for the mass de- 
termination. For the eclipsing X-ray source 4U 1700-37 with 
the 06.5 laf+ companion HD 153919 (Jones et al. 1973; Mason 
et al. 1976) no X-ray pulsations have been detected, although 
the compact object most likely is a neutron star (Reynolds et al. 
1999; Van der Meer et al. 2005). The absence of X-ray pulsations 
prohibits the accurate determination of the orbital parameters of 
the neutron star, and thus its mass. 

Van Kerkwijk et al. (1995a) present an analysis of the 
neutron-star mass determinations for these systems hosting an 
X-ray pulsar and conclude that the accuracy of the (then) avail- 
able observations does not allow to discriminate between one 
"canonical" neutron-star mass or a mass distribution. Recent 
analyses of the radial-velocity curve of the wind-fed system 
Vela X-1 (Barziv et al. 2001; QuaintreU et al. 2003) with its 
B0.5 lb companion (Hiltner et al. 1972; Vidal et al. 1973) 
have shown that the neutron star in this system has a mass of 
1.86 + 0.16 M©. Such a high neutron-star mass provides an im- 
portant constraint on the EOS at supra-nuclear density. 

Since the work of Reynolds et al. (1992), included in 
the analysis by Van Kerkwijk et al. (1995a), no new optical 
spectroscopy of the BO supergiant companion (QV Nor) of 
4U 1538-52 has been reported in literature. Van Kerkwijk et al. 
(1995a) Ust 1.06:^°;^i M© for the mass of 4U 1538-52. 

We have obtained VLT/UVES spectra of the three Roche- 
lobe overflow systems SMC X-1, LMC X-4 and Cen X-3. The 
orbital parameters (Table 1) of their X-ray pulsars are accurately 
known, based on X-ray pulse time delay measurements. The X- 
ray pulsars in these systems have short spin periods (seconds) 
compared to those in wind-fed systems (minutes), as the mass- 
and angular-momentum accretion rate in Roche-lobe overflow 
systems is much higher than in wind-fed systems. The photo- 
metric light curves indicate that in all three systems an accre- 
tion disc is present (Tjemkes et al. 1986; Heemskerk & Van 
Paradijs 1989). The X-ray eclipse duration is best measured in 
hard X-rays, since at lower energies the eclipses are systemat- 
ically longer due to soft X-ray absorption by the stellar wind 
of the OB companion. The eclipse duration can thus be used to 
determine the radius of the OB companion. 
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Table 2. Observing log of the three observed systems. For each observation the Modified Julian Date (MJD) corresponding to the mid-exposure 
time is given, the orbital phase (0) and the signal-to-noise ratio (S/N) per resolution element of the three different CCDs centered at 4270 A for 
the blue CCD, 5150 A for the redl CCD and 6100 A for the red2 CCD. 



SMC X-1 


LMC X-4 


Cen X-3 


MJD 


orbital 




S/N 




MJD 


orbital 




S/N 




MJD 


orbital 




S/N 




(days) 


phase(0) 


blue redl red2 


(days) 


phase(0) 


blue redl red2 


(days) 


phase(0) 


blue redl red2 


52187.154 


0.537 


58 


75 


69 


52214.187 


0.388 


54 


69 


58 


52271.300 


0.462 


33 


67 


77 


52214.090 


0.458 


60 


80 


71 


52224.249 


0.533 


43 


58 


50 


52287.229 


0.095 


22 


50 


64 


52224.224 


0.062 


49 


69 


62 


52225.255 


0.247 


49 


66 


57 


52292.273 


0.512 


27 


57 


68 


52225.211 


0.315 


58 


78 


68 


52226.212 


0.927 


48 


67 


56 


52298.169 


0.337 


32 


63 


75 


52242.059 


0.644 


41 


61 


54 


52237.314 


0.809 


54 


72 


63 


52309.323 


0.681 


35 


76 


81 


52243.087 


0.908 


66 


77 


70 


52242.282 


0.337 


55 


71 


61 


52313.150 


0.515 


37 


74 


78 


52244.168 


0.186 


67 


81 


74 


52243.113 


0.927 


54 


62 


55 


52321.163 


0.354 


29 


64 


72 


52247.185 


0.961 


43 


58 


51 


52244.193 


0.694 


59 


67 


55 


52322.158 


0.831 


31 


65 


73 


52256.112 


0.255 


53 


66 


59 


52245.213 


0.418 


67 


73 


65 


52326.186 


0.761 


38 


77 


82 


52258.094 


0.764 


58 


69 


62 


52246.085 


0.037 


44 


53 


45 


52327.277 


0.284 


37 


77 


85 


52270.052 


0.837 


61 


74 


68 


52256.086 


0.138 


32 


40 


35 


52328.241 


0.746 


41 


82 


88 


52271.059 


0.095 


52 


65 


60 


52259.094 


0.274 


49 


59 


50 


52329.260 


0.234 


38 


76 


86 


52285.112 


0.706 


47 


61 


56 


52261.163 


0.743 


58 


71 


61 













2.1. SMC X-1 

The BO supergiant Sk 160 (V - 13.3 mag) is the companion 
to the eclipsing X-ray pulsar SMC X-1 (Schreier et al. 1972a; 
Liller 1973), located in the "wing" of the Small Magellanic 
Cloud at a distance of 60.6 kpc (Hilditch et al. 2005). The spin 
period of the pulsar is 0.71 s and the orbital period is 3.89 d, 
which is decaying on a timescale of 3 x 10^ yr due to tidal inter- 
action (Levine et al. 1993). A super-orbital, though not strictly 
periodic variation of ~ 60 d is present in the system, most likely 
due to a precessing tilted accretion disc (Wojdowski et al. 1998; 
Clarkson et al. 2003). 

The most recent determination of the radial-velocity orbit 
of Sk 160 has been performed by Val Baker et al. (2005). 
Optical spectra covering the wavelength range 4300-5100 A 
were obtained with the grating spectrograph on the 1.9 metre 
Radcliff telescope at the Sutherland Observatory, with a resolv- 
ing power R ~ 4000. The majority of the 56 usable spectra were 
secured during one week of observations in September 2000. 
Based on a cross-correlation analysis similar to the one used 
by Reynolds et al. (1993), a radial-velocity amplitude of K^pt = 
18.0±1.8 km s~^ was measured, which becomes 21. 8+1. 8 km s~^ 
when taking the effects of X-ray heating into account; the rest- 
frame (y) velocity is 174 km s ' . To simulate the effects of X-ray 
heating, a model is used that generates velocity corrections based 
on contributions from different elements of the projected stellar 
disc. The models do not take into account the presence of an 
accretion disc, which may well reduce the effect of X-ray heat- 
ing (Van Kerkwijk et al. 1995a). The results of Val Baker et al. 
(2005) are significantly different from the results obtained by 
Reynolds et al. (1993) who arrive at K^pt = 27.5 + 1.9 km s"' fol- 
lowing a similar procedure. According to Val Baker et al. (2005), 
this discrepancy could be due to the limited phase coverage of 
the dataset of Reynolds et al. (1993) and to the fact that Reynolds 
et al. (1993) assume a significantly higher value for Lx when de- 
termining the non-Keplerian corrections. The latter would, how- 
ever, not explain the difference in Kgpt before applying the X-ray 
heating corrections. Val Baker et al. (2005) derive lower (edge- 
on system) and upper limits (Roche-lobe filling system) to the 
mass of SMC X-1 of 0.91±0.08 M© and 1.21±O.1OM0, respec- 
tively. The mass of the optical companion is around 16.7 M© in 
both cases. 



2.2. LMC X-4 

After the first detection of LMC X-4 by the Uhuru satellite 
(Giacconi et al. 1972), the binary nature of its optical coun- 
terpart was confirmed by Chevalier & llovaisky (1977). The 
V = 14.0 mag 08 111 companion (Sanduleak & Philip 1976; 
Kaper et al., to be submitted) is in a 1.41 d orbit (Li et al. 1978; 
White 1978), which is decaying on a timescale of ~ 500 000 yr 
(Levine et al. 2000). The optical light curve shows ellipsoidal 
variations and a super-orbital period of ~ 30 d due to a precess- 
ing accretion disc (Heemskerk & Van Paradijs 1989). The X-ray 
light curve includes regular eclipses as well as a pronounced flux 
modulation of a factor ~ 60 with a period of 30.5 d (Lang et al. 
1981). This long-term variation is attributed to the precessing 
accretion disc. Kelley et al. (1983) discovered the 13.5 s X-ray 
pulsations of LMC X-4. 

Chevalier & llovaisky (1977) reported on photographic spec- 
tra obtained with the 1.5m ESO telescope from which they 
derived radial-velocity variations with an amplitude of 475 ± 
25 km s~' for the Hell 4686 A line, and a phase dependence 
suggesting an origin near the X-ray source. Hutchings et al. 
(1978) collected 18 spectrograms using the Cassegrain image- 
tube spectrograph of the CTIO 4m telescope in November 1977, 
with an effective spectral resolution ofR ~ 3000. For the hydro- 
gen fines (HjS to H9 , and an empirically determined correction 
to the blended line) they derive a radial-velocity amplitude 
of A'opt =50 + 5 km s"' ; for the He I lines Ko^t = 60 ± 9 km s ' is 
measured. For the He ii 4686 A line Hutchings et al. (1978) de- 
rive Kopt = 498 ± 14 km s~\ with a phase difference compared 
to the H and He i absorption lines of 0.79 f orb, consistent with 
ChevaHer & llovaisky (1977). Kelley et al. (1983) combined 
the radial-velocity data of the hydrogen absorption lines pre- 
sented by Hutchings et al. (1978) with measurements by Petro 
& Hiltner (1982) and arrive at Kopt = 37.9 ± 2.4 km s-^. Van 
Kerkwijk et al. (1995a) use Kopt = 38 + 5 km s~^ and obtain 
Mx = 1.47:0;^ Mo and Mopt = 15.8:|3 m©. 

2.3. Cen X-3 

Cen X-3 was discovered by Chodil et al. (1967) and became the 
first detected binary X-ray pulsar (Giacconi et al. 1971; Schreier 
et al. 1972b). The V = 13.3 mag optical counterpart V779 Cen 
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Fig. 1. Normalised spectra in the wavelength range 3580-4500 A of SMC X-1 at orbital phase = 0.06 (top spectrum), LMC X-4 at orbital 
phase (p = 0.04 (middle spectrum) and Cen X-3 at orbital phase = 0.09 (bottom spectrum), respectively. The line identifications are shown 
above the spectrum of SMC X-1. Note that the reddening of Cen X-3 (E{B -V) ~ 1.4) clearly affects the S/N in this wavelength range. 
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Fig. 2. Normalised spectra in the wavelength range 4625-5585 A of SMC X-1 at orbital phase <f> = 0.06 (top spectrum), LMC X-4 at orbital 
phase (p = 0.04 (middle spectrum) and Cen X-3 at orbital phase (p = 0.09 (bottom spectrum), respectively. The line identifications are shown 
above the spectrum of SMC X-1. 
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above the spectrum of SMC X-1. 
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Fig. 4. Example of gaussian fits for the Hei line at 4387.93 A (upper 
panel), the Hei line at 4713.17 A (middle panel) and the Hei line at 
5047.74 A (lower panel) in the spectrum of the B supergiant companion 
to SMC X-1. The spectra and the gaussian fits are shown at four dif- 
ferent quadratures. To better visualise the fits the spectra are smoothed. 
Note that the gaussian profile fits represent the rotationally broadened 
line profiles very well. Furthermore, the orbital motion of the B super- 
giant companion is easily observed. 



was identified by Krzeminski (1974), an 06-7 II-III star (Ash 
et al. 1999) in a 2.09 d circular orbit with the 4.84 s X-ray pul- 
sar. The optical light curve indicates the likely presence of an 
accretion disc, but no strong evidence is found for X-ray heating 
(Tjemkes et al. 1986). The X-ray light curve includes episodes 
of high and low X-ray flux with a characteristic timescale of 
120-165 d (Priedhorsky & Terrell 1983; Paul et al. 2005). 

Based on photographic spectra, Hutchings et al. (1979) 
derive Kopt = 24 ± 6 km s"\ confirmed by Aslanov & 
Cherepashchuk (1982) who also report Aopt = 24 ± 6 km s"'. 
The most recent radial- velocity measurements of V779 Cen are 
presented by Ash et al. (1999) who determine two very dif- 
ferent values of /Topt based on two datasets obtained with the 
4m Anglo-Australian Telescope and the RGO spectrograph. The 
wavelength range of these spectra is 4300^700 A; the spectral 
resolution R ~ 3000. Ash et al. (1999) discard the results of flie 
first dataset and arrive at A'opt = 24.4 + 4. 1 km s"' . The resulting 
neutron-star mass is Mx - 1.21 + 0.21 M© and the mass of the 
O-type companion Mopt - 20.5 + 0.7 Mq. 
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Fig. 5. Example of radial-velocity curves obtained from the HI 3 line at 
3734.37 A. Some datapoints are shown twice to better visualise trends 
with orbital phase. The error bars indicate 1 cr errors. The lower panels 
show the residuals of the fit. 



3. Observations 

We have obtained high-resolution (R ~ 40 000) spectra of the 
three systems with UVES (Dekker et al. 2000) on the VLT in ser- 
vice mode in the period October 2001 to March 2002 at Paranal, 
Chile. The total exposure time was 17.6 h spread over 13 ex- 
posures of 1400 s of SMC X-1, 13 exposures of 2000 s of 
LMC X-4 and 12 exposures of 1600 s of Cen X-3. The in- 
strument was used with standard setting "390-1-564" and a slit 
width of 1.0". This yields a wavelength range of 3580^500 A 
in the blue arm and 4625-5585 A and 5680-6640 A in the red 
arm. To determine the orbital phase of the systems we used the 
ephemeris of Wojdowski et al. (1998), Levine et al. (2000) and 
Nagase et al. (1992) for SMC X-1, LMC X-4 and Cen X-3, 
respectively (see Table 1). The log of observations is listed in 
Table 2. 
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Fig. 6. The upper panels show the radial- velocity amplitude (A'opt in km s"') measured for a given line as a function of its mean line equivalent 
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the (mean) radial-velocity amplitude. These selections show less spread in A^opt- 



In order to reduce the data to normalised spectra we used the 
UVES pipeline (version 1.2.0) and the ESO reduction package 
MIDAS (version 03SEPpIl.l). All raw echelle frames were bias 
and flatfield corrected. Subsequently the different orders were 
extracted using the optimal extraction routine available within 
the UVES pipeline. High signal-to-noise (S /N ^ 50) spectra 
extracted with this routine show a ripple effect. In our spectra this 
effect is only marginally present in some spectra of SMC X-1 
in the wavelength range 5100-5585 A. This range contains the 
He II 541 1 .53 A line, but we do not use it for the determination 
of the radial-velocity amplitude of the system (Sect. 4.2). 

Cosmic ray hits were removed by rejecting the affected 
wavelength bins and subsequent interpolation. After this, the 



spectra were normalised by fitting the continuum with a spline 
over the whole wavelength range of one spectral arm. The nor- 
malised spectra of the three systems are shown in Figs. 1 to 3. 



We verified the long-term stability of UVES by measuring 
the position of the interstellar lines of Can K at 3933.66 A, 
Can H at 3968.47 A, Nai Dl at 5895.92 A and Nai D2 at 
5889.95 A for all spectra. This resulted in a deviation of less 
than 1 km s ' throughout the whole observing period of each 
system, i.e. several months. 
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Table 3. List of lines that have been fitted with a gaussian. The rest wavelength of the line is given as well as the mean equivalent width (EW, in 
Angstrom) with its Icr error. The magnitude of the variations in line EW is expressed as the ratio of the standard deviation of the line EW to the 
error on the mean EW. Rejected lines are indicated with a remark: lines marked with an "B" are blends; a "V" stands for lines exhibiting strong 
EW variations (77 > 9); a number (e.g. 0.5) indicates a deviation in the EW at that orbital phase. Electronic tables that contain the radial-velocity 
amplitude, EW and IWHM variations for each line are available as on-line material at the CDS. 
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4. Spectral Analysis 

To obtain a radial-velocity measurement often the complete 
spectrum is cross-correlated with a template spectrum. This ap- 
proach has many advantages when using spectra with relatively 
low spectral resolution and poor signal-to-noise. In our case the 
spectra are of such high quality that the radial-velocity ampli- 
tude can be determined for each line separately. The advantage 
of such a strategy is that it is possible to assess the influence of 
possible distortions due to e.g. X-ray heating and gravity darken- 
ing, as in these systems the OB star is irradiated by a powerful X- 
ray source (Lopt ^ Lx) and is filling its Roche-lobe. Furthermore, 
the extended OB-star wind is focused into a shadow wind which 
possibly produces a strong shock (a so-called photo-ionisation 
wake, see e.g. Blondin 1994; Kaper et al. 1994) where the fast 
shadow wind catches up with the stagnant flow inside the X-ray 
ionization zone. The shadow wind and photo-ionization wake 
induce orbital modulations of spectral lines formed in the stel- 



lar wind (i.e. strong spectral lines such as the first lines of the 
Balmer series and the strongest heUum lines). 

Figures 1 to 3 show that the spectra contain mostly lines 
that are identified with transitions from H and He i. Only a few 
Hen fines and some metal fines are detected, consistent with the 
modest metallicity of the Magellanic Clouds and OB supergiant 
spectral types. We show the spectra observed near X-ray eclipse 
(orbital phase <p ~ 0.0). Some of the lines are blended or show 
a sfight asymmetry. A comparison of the spectra obtained at dif- 
ferent orbital phase reveals that several lines also vary in line 
strength. Still, Fig. 4 demonstrates that many lines are well rep- 
resented by a gaussian profile (i.e. as one would expect for a 
rotationally broadened profile). Therefore, gaussian profiles are 
fit to all individual lines to determine the radial-velocity curve 
(Sect. 4.1). Subsequently, the observed line profiles are exam- 
ined on asymmetry and variations with orbital phase (Sect. 4.2). 



10 



A. van der Meer et al.: Determination of the mass of the neutron star in SMC X-1, LMC X-4 and Cen X-3 



0.6 

g 0.4 
S 

LlI 

0.2 



He II 541 1.53 A 
7) = 42.5 I 



0.8 

S 

Ld 

0.6 

0.4 
2 

S 1.8 
S 

Ld 

1.6 

1.4 

1 

g 0.8 
S 

LlJ 

0.6 
0.4 



He I 4471.50 A 
7j = 12.3 



— 1 1 r- 



i 

I I 



i- 



He 3970.07 A 
7) = 8.5 



— 1 1 r- 



H12 


3750. 


1 1 i i 1 1 1 i i i 1 1 

5 A 


V = 


6.6 




h 


» i 

r 




1 


1 1 






0.5 1 

Orbital Phase (0) 



1.5 



Fig. 7. Equivalent width (EW) variations with orbital phase (^) ob- 
served in SMC X-1. The error bars indicate 1 cr errors and the dotted 
Une the average value. Some datapoints are shown twice to better visu- 
alise trends with orbital phase. The value of the defined line variability 
parameter ?] is indicated for each line. The upper panel shows the large 
variations in the Hen line at 5411.53 A, the second panel shows the in- 
verse behaviour of the He i line at 4471.50 A, the third panel shows the 
variations of He at 3970.07 A. A stable line is shown in the lower panel, 
i.e. the H12 line at 3750.15 A. The EW ratio of the Hen and He i lines 
is sensitive to Teff and can be used for spectral classification; apparently, 
the spectral type of SMC X-1 varies with orbital phase. The variations 
are consistent with being caused by X-ray heating of the stellar surface 
facing the X-ray source. 



4.1. Radial-velocity curves 

We determine the Une centre, and thus the Doppler shift with 
respect to the heliocentric restframe, by fitting the profile with 
a gaussian. The gaussian sets the full-width at half maximum 
(FWHM), the central hne depth and the central wavelength of 
the profile, i.e. three free parameters. A minimahzation pro- 
cedure delivers the best fit gaussian profile and defines the ac- 
curacy of the fit parameters. Emission-line profiles, such as Ha, 
are not included in the fitting procedure. 

The orbital parameters are accurately known from X-ray 
pulse time delay measurements (see Table 1). We assume the or- 
bit to be circular, since in all cases the eccentricity e < 0.008. 
The remaining free parameters describing the radial-velocity 
curve are the radial-velocity amplitude, Ko^t, and the restframe 
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at 4340.46 A (upper panel) and the Hei hne at 4387.93 AA (lower 
panel) of SMC X-1 as a function of orbital phase (0) in km s~^ The 
error bars indicate 1 cr errors and the dotted line the average value. 
Some datapoints are repeated once to better visualise trends with or- 
bital phase. A maximum in FWHM is clearly seen around (p ~ 0.25; a 
second, though much less pronounced maximum is visible at ^ ~ 0.75. 
Such FWHM variations are also detected in other lines, as well as in the 
systems LMC X-4 and Cen X-3. 



velocity of the system, Vy. In principle, a small shift in orbital 
phase could be present due to the inaccuracy of the orbital pe- 
riod, the period derivative, and mid-eclipse time (Table 1). The 
measured orbital phase shifts are not significantly different from 
zero, but are slightly larger than the phase shifts one may expect 
based on the accuracy of the ephemeres of these systems. The 
phase shift should be the same for all lines and has to be fixed to 
the average value during a second iteration. 

A radial- velocity curve is obtained for each individual hne; 
Fig. 5 displays for each system a radial-velocity curve represen- 
tative for the spectral lines used to measure the radial-velocity 
amphtude. Note that the data points were obtained from several 
orbits of the system; the dozen spectra per system are evenly 
distributed with orbital phase, thanks to the service-mode obser- 
vations allowing to obtain spectra spread over a period of more 
than one month. 

Fig. 6 indicates that A'opt and Vy show quite some dispersion 
when comparing the radial-velocity curves of individual lines. 
Kofx shows a dependence on line strength and line variability, 
diagnostics that we will use to reject lines when determining the 
mean radial- velocity amphtude used to calculate the mass of the 
neutron star. 



4.2. Line Selection 

As our observational strategy is aimed at the derivation of the 
radial-velocity curve based on individual lines, we select only 
lines that are well identified. Also, all lines that are (partially) 
blended, are rejected. In Table 3 this is indicated with a "B" as a 
remark. 

In these high-mass X-ray binaries it is expected that at least 
some lines in the OB-star spectrum are affected by the presence 
of the X-ray pulsar companion (see e.g. Van Paradijs et al. 1978; 
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Table 4. Final selection of lines that is used for the determination of the radial-velocity amplitude (A'opt) for each system. The three fit parameters 
are listed, i.e. the system velocity (r,,,). A",,,,,, and a phase shift (A(fi). The goodness of the fit is expressed asx^-- All errors are 1 a. 
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Reynolds et al. 1993). In principle, one could model the pho- 
tospheric line profiles of the OB supergiant companion by cal- 
culating a grid of many (thousands) surface elements, adopting 
for each surface element an intrinsic line profile, and obtain- 
ing the integrated line profile by a weighted integration of the 
(mainly by stellar rotation) Doppler-shifted intrinsic profiles of 
all visible surface elements. Such computationally demanding 
techniques have been successfully explored in stellar pulsation 
studies (e.g. Schrijvers et al. 1997) and in modeling photometric 
Ughtcurves of HMXBs (e.g. Heemskerk & Van Paradijs 1989). 
Also, variations in the intrinsic line profiles due to variations of 
surface temperature and gravity have been included (Schrijvers 
& Telting 1999), so that one could apply such a method to the 
case of an OB supergiant irradiated and deformed by a close and 
compact X-ray source. Although we have now made an advance 
in analysing spectra of the OB companions to X-ray pulsars by 
studying individual lines rather than cross-correlating complete 
spectra, we consider this modeling effort beyond the scope of 
the present paper, even though it has the potential to deliver very 
important information on the interpretation and analysis of the 
obtained spectra. 

A first step in this direction has been undertaken by 
Abubekerov et al. (2004); they model the line formation process 
taking X-ray heating and gravitational darkening into account. 
As the luminosity distribution and shape of the star are altered, 
the measured radial velocity, as well as the line equivalent width, 
depend on the depth of the line-forming region and orbital phase. 
Abubekerov et al. (2004) show that these effects will result in 
a reduction of the derived radial- velocity amplitude by several 
km s 

That these effects manifest themselves in our observations 
is nicely demonstrated by the large EW variations of the He n 
line at 5411.53 A of SMC X-1 (see Fig. 7). This line varies 
with about a factor ~ 5 in EW and reaches a maximum EW near 
(p ~ 0.5. The He: line at 4471.50 A (and also Hei 5875.66 A) 
shows the opposite behaviour and has a minimum strength at this 
orbital phase (see Fig. 7). Since a He i / Hen line ratio is sensi- 
tive to Tgff it can be used for spectral classification; however, 
traditionally the ratio of Hex 4471 over Heii 4541, accessible 



in the blue spectrum, is used for spectral classification (Conti 
& Alschuler 1971; Lennon et al. 1993; Mokiem et al. 2005). 
Line ratios involving the Hei 5876 and Hen 5412 lines have not 
(yet) been caUbrated (Mokiem, priv. comm.). The observed vari- 
ations in line ratio indicate a higher T^g at the side of the OB 
star facing the X-ray source, and is evidence for X-ray heating 
in SMC X- 1 . As a consequence, the spectral type of the optical 
companion to SMC X- 1 changes as function of binary aspect 
angle. A similar effect has been observed in the optical spectrum 
of V779 Cen, the 06 companion of Cen X-3 (Hutchings et al. 
1979). 

The derived value of the radial-velocity amplitude shows a 
dependence on line strength (Fig. 6). Furthermore, the observed 
variations in EW also increase with line strength. As illustrated 
above, these variations likely reflect distortions of the line form- 
ing region due to e.g. X-ray heating and the (disturbed) stel- 
lar wind. Since K^pt should be a unique value, we investigate 
whether a selection criterion can be defined to reject lines that 
are affected by these distortions and thus do not yield a sound 
measurement of ^opt- 

To measure these distortions we apply a velocity moment 
analysis, with which one can determine the equivalent width 
(EW), central velocity, standard deviation (cr) and skewness (r) 
of a given line. The n* moment (//„) of a distribution /(v) in 
velocity v is given by: 



/ V" /(v) dv 
Jmdv ' 



and 



= j f{v)dv. 



(1) 



(2) 



For such a distribution the EW is proportional to /zq and the cen- 
tral velocity x.o The standard deviation {&) and skewness (r) 
are often defined as: 



/ iv-fxifmdv 



(3) 
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Table 5. Weighted mean values of the radial velocity for each spectrum, indicated hy its orbital phase ((p) for SMC X-1, LMC X-4, and Cen X-3. 
All errors are 1 cr. 



SMC X-1 


LMC X-4 


Cen X-3 





< Vrad > 








< Vrad > 




(kms-') 




(kms-') 




(km s-') 


0.062 


8.42 ± 1.06 


0.037 


6.45 ± 1.50 


0.095 


14.96 ± 2.97 


0.0% 


19.04 ± 1.02 


0.138 


29.28 ± 2.23 


0.234 


22.87 ±2.10 


0.186 


16.58 + 0.77 


0.247 


40.06 ± 1.38 


0.284 


31.47 ± 1.98 


0.255 


24.42 ± 1.17 


0.274 


30.03 ± 1.60 


0.337 


14.12 ±2.59 


0.315 


14.70 ± 0.95 


0.337 


24.84 ± 1.24 


0.354 


18.43 ±2.75 


0.458 


4.98 ± 0.84 


0.388 


25.04+ 1.33 


0.462 


-10.75 + 2.86 


0.537 


-9.37 + 1.01 


0.418 


23.23 + 1.10 


0.512 


-21.61 + 2.42 


0.645 


-18.46 ± 1.24 


0.533 


-14.69+1.54 


0.681 


-26.12 ±2.19 


0.707 


-19.07 ± 1.08 


0.694 


-30.90 ± 1.21 


0.746 


-18.81 ± 1.76 


0.765 


-16.47 ±0.89 


0.743 


-34.06 ± 1.17 


0.761 


-23.29 ± 2.07 


0.837 


-19.83 ±0.87 


0.809 


-30.86 ± 1.16 


0.831 


-24.32 ± 2.34 


0.909 


-8.22 ± 0.77 


0.927 


-17.96 ± 1.19 


0.515 


-14.55 ± 1.98 


0.961 


-1.85 ± 1.09 


0.927 


-17.95 ± 1.12 







and 

in which the skewness (r) is a measure of the asymmetry of the 
distribution. For a gaussian profile the skewness t is zero; if t 
the spectral line is not well represented by a gaussian. If so, 
one can not use a gaussian to fit the Une in order to measure 
its radial velocity, as we did in Sect. 4. 1 . It turns out that for the 
lines that are not blended, r is consistent with being equal to zero 
within the error. Note, however, that the error on t is too large to 
measure any significant deviations from zero, because it depends 
on the error on fi and a to the third power (see Eq. 4). The EW 
and central velocity (first moment) are much more accurately 
determined and are consistent with the respective values derived 
from the gaussian profile fits. 

For each line we determine whether the hue EW varies sig- 
nificantly by comparing the deviation in EW to the error on the 
mean EW (Table 3). Note that these values are not equal; the 
error on the mean depends on the error in EW of the individual 
spectra, while the standard deviation is the spread in the distribu- 
tion of EWs. Since the strongest lines are formed in the outer lay- 
ers of the stellar photosphere and/or in the extended stellar wind 
these lines are expected to be most affected by X-ray heating, 
gravity darkening, etc., and will thus show intrinsic variations 
when the system revolves. The observed trend in Ko^t with Vy 
(Fig. 6) also indicates that the stronger lines are formed further 
out in the stellar atmosphere and wind (an effect called Balmer 
progression, see e.g. Crampton et al. 1985; Abubekerov et al. 
2004). We define a fine variabihty parameter rj to formulate a 
selection criterion (Table 3). This parameter is defined as the ra- 
tio of the standard deviation of the EW variations to the error on 
the mean EW. The main motivation behind this definition is that 
the EW is in principle not sensitive to radial-velocity variations, 
i.e. the key parameter that we want to measure in the spectra. 
However, t] turns out to be an accurate probe of intrinsic Une 
profile variability. Other methods to measure line profile vari- 
ability (e.g. the temporal variance spectrum analysis method in- 
troduced by FuUerton et al. 1996) are sensitive to radial- velocity 
variations. 

We select 77 = 9 as the value above which lines are rejected. 

These lines are indicated with a "V" in Table 3. The threshold 
value for 77 is chosen arbitrarily, but is a reproducable and ob- 



jective means to quantify line profile variations. As it is an av- 
eraged quantity, the rj threshold does not reject lines exhibiting 
only modest EW variations and having relatively large errors on 
the measured EWs (e.g. the He at 3970.07 A shown in Fig. 7). 
On the other hand, it does reject some lines that have highly ac- 
curately determined EWs and that exhibit hardly any EW varia- 
tions with orbital phase. To test the impact of the chosen value 
for the T] threshold on the obtained value of the radial-velocity 
amplitude, we have evaluated the results for different values of 
77. Including all lines for which ?/ < 8 then Kopt is 20.8 ± 1.5, 
35.7 ± 1.8, and 27.5 ± 2.3 km s"^ for SMC X-1, LMC X-4, and 
Cen X-3, respectively. Similarly, we obtain 20.0± 1 .2, 34.2± 1 .3, 
and 26.2 + 2.7, respectively, if we select the lines with 77 < 13. As 
expected, we obtain a lower value of Kgpi for a higher value of 77, 
and vice versa, while the error on the result increases when in- 
cluding lines that show more intrinsic variations. Still, the values 
for A'opt agree within the error for the applied range in 77. 

We now fix the threshold to 77 = 9; to ensure that the lines 
with 77 < 9 and small variations concentrated around (p ~ 0.0 or 
4> ~ 0.5, are rejected, we mark these in Table 3 hsting the orbital 
phase at which these variations are concentrated. Note that most 
of these lines are already marked with a "V", i.e. rejected based 
on the 77 threshold. 

The full width at half maximum (FWHM) of the line also 
varies as a function of orbital phase 0. This behaviour is espe- 
cially visible in the Balmer series lines of hydrogen, most promi- 
nent in the stronger lines. Figure 8 shows that the FWHM in- 
creases when the system revolves from X-ray eclipse to <f> ~ 
0.25, where it reaches a maximum before declining again. At 
(p ~ 0.75 another, though smaller increase in FWHM is de- 
tected. The helium (and some other metallic) lines also show 
FWHM variations, but less pronounced and not as periodic as 
observed in the hydrogen Unes. Furthermore, these FWHM vari- 
ations are best seen in SMC X-1 and LMC X-4, but are less 
clear in Cen X-3. 

Apparently, if one would use these lines to determine the 
(projected) stellar rotation velocity, one would arrive at a larger 
value of V sin ; when the system is looked upon from a side view. 
This may be due to the elongated shape of the star as it is filling 
its Roche lobe, as evidenced by the observed ellipsoidal varia- 
tions. The dependence on line strength would be explained by 
the fact that the line forming region is further out in the atmo- 
sphere when the line is stronger. This would, however, not ex- 
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plain the difference in amplitude of this effect observed between 
-0.25 and -0.75. 

Differences in spectral appearance when comparing spectra 
obtained at - 0.25 and - 0.75 are well known to occur 
in spectroscopic binaries. The "Struve-Sahade" effect (Struve 
1937; Sahade 1962) is the apparent strengthening of the sec- 
ondary spectrum of a hot binary when the secondary is approach- 
ing and the corresponding weakening of the lines when it is 
receding (see Gies et al. 1997 for an observational overview). 
The cause of this effect may be the presence of a gas stream 
(bow shock, wind collision) trailing the secondary in its orbit 
(Sahade 1959; Gies et al. 1997). Hydrodynamical simulations of 
SMC X-1 by Blondin (1994) indicate that a collision of wind 
material from the shadow wind with material in the X-ray ion- 
isation zone is present in the system. Perhaps that this shocked 
material introduces the difference in amplitude observed in the 
FWHM of the strong Balmer lines that are formed in the stellar 
wind. 

As we do not have a clear explanation for these variations 
in FWHM, we investigate the possible influence of this effect on 
the determination of Kopi- Note that the stronger Balmer lines are 
already excluded from the radial-velocity analysis on the basis of 
the EW variations. We fix the FWHM on its mean value, its max- 
imum and its minimum and refit the lines. It turns out that the 
measurement of the centre of the line profile is not affected much 
by these FWHM variations. The derived radial- velocity ampli- 
tude, Kopt, is the same within its errors in all cases. Therefore, 
we decide not to exclude more lines based on a FWHM variation 
criterion (most of the lines showing this effect were excluded on 
other grounds anyway). 



4.3. Radial-velocity amplitude 

Table 4 shows the final selection of Unes for which good fits to 
the radial-velocity curve are obtained, resulting in a measure- 
ment of Kopt and Vy. Since the phase shift should be equal for all 
lines in one system, we refit all lines with the phase shift fixed 
to the weighted average, i.e. A(f> is 0.026, -0.003, and 0.065 for 
SMC X-1, LMC X-4, and Cen X-3, respectively. Fixing fliese 
values does not influence the other parameters much; they are 
the same within their errors. 

To avoid sytematic errors while determining the final mean 
value of ^Topt we shift all determined line centres to equal Vy. The 
weighted mean values of Vy are 191+6 km s"', 306 + 10 km s"' 
and 32 ± 13 km s ^ for SMC X-1, LMC X-4 and Cen X-3, 
respectively. Then we calculate the weighted mean of the radial 
velocity for each spectrum (see Table 5). The average dataset we 
fit in the same way as the individual lines, which results in the 
three radial-velocity curves shown in Fig. 9. The goodness of the 
fits with respect to the number of degrees of freedom (d.o.f.) are 
xj/d.o.f = 9.S/l2,x^r/d-0-f = 9.1/12 and ;f?/d.o.f = 7.2/11 for 
SMC X-1, LMC X-4 and Cen X-3, respectively. The error bars 
indicate 1 cr errors multiplied hyxr to obtain ax^ - 1 . The resid- 
uals to the fit do not show any further evidence for systematic 
effects with orbital phase, as is the case for e.g. Vela X-1 (Van 
Kerkwijk et al. 1995b; Barziv et al. 2001). The final values of 
Kopt are 20.2 + 1 .1 km s"', 35.1 ± 1.5 km s"' and 27.5+2.3 km s"' 
for SMC X-1, LMC X-4 and Cen X-3, respectively. The accu- 
racy of the determination of Kgpt has been significantly improved 
(by a factor 2-4) compared to previous measurements (and com- 
parable to and consistent with Val Baker et al. (2005) in the case 
of SMC X-1). These values are subsequently used to calculate 
the mass of the optical companion and the X-ray source. 




10 

-10 



11 
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Orbital Phase {</>) 



1.5 



Fig. 9. Radial- velocity curves obtained from weighing the radial veloc- 
ities of all selected lines per observation. Some datapoints are shown 
twice to better visuaUse trends with orbital phase. The error bars indi- 
cate 1 cr errors multiplied by Xr to obtain ax^ = 1. Note that Vy has 
been set to zero. 

5. The neutron star masses 

In order to measure the mass of the neutron star and its optical 
companion we apply the mass function. For an orbit with eccen- 
tricity e it can be shown that this is defined as: 



Mopt 
and 



2nG sin^ i 



(1 + qf 



27rGsin i 



(5) 



(6) 



where Mopt and Mx are the masses of the optical component 
and the X-ray source, respectively, Kopt and Kx are the semi- 
amphtude of the radial- velocity curve, P is the period of the orbit 
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Fig. 10. Selection of lines for which we measured the rotational broadening. Each diagram includes a solid hne representing the model that we 
appUed; the dotted line represents the continuum. To better visualise the model, the spectra are smoothed with 7 resolution elements. The top 
panels show the modelled spectral lines of SMC X- 1 , the second and third row include lines of LMC X-4 and the two bottom rows show lines of 
Cen X-3. 



and i is the inclination of the orbital plane to the line of sight. The inclination of the system we follow the approach of Rappaport 
mass ratio q is defined as: & Joss (1983), who showed that: 



Opt 



^Opt 



(7) 



smi 



COS d« 



(8) 



The values for Kx and P can be obtained very accurately where /?l is the Roche-lobe radius of the optical component, jS 

from X-ray pulse timing delay measurements (Wojdowski et al. is the ratio of the radius of the optical component to /?l (i-e. a 

1998; Levine et al. 2000; Nagase et al. 1992). The VLT/UVES Roche-lobe filling factor), a is the separation of the centres of 

spectra provide a value for Kopt. For the determination of the mass of the two components, and 0e is the semi-eclipse angle 



A. van der Meer et al.: Determination of the mass of the neutron star in SMC X-1, LMC X-4 and Cen X-3 



15 





1.5 0.6 0.8 1 1.2 1.4 0.4 0.6 0.8 1 

Q Q 



Cen X-3 


1 1 1 1 1 1 1 1 1 


, , , 1 





60 70 ^ 80 60 65 70 ^ 75 80 60 70 80 90 

Inclination C) InclInatlonC) Inclination Q 





1 1 1 1 1 
Cen X-3 


1 1 1 1 1 


--^-^ , . , 


, ■ 1 . 



0.55 0.6 0.65 




12 14 16 18 20 



14 ie 

Wo (Mo) 



18 15 



20 
Mo (Ms) 




1 1.2 




'I'll 
Cen X-3 


1 1 1 1 1 







1.2 1.4 



Fig. 11. Resulting probability distributions from the Monte-Carlo simulations for SMC X-1 (left column), LMC X-4 (middle column) and 
Cen X-3 (rigth column). Shown are the ratio of the rotational and orbital frequency (£1), the incUnation angle (?'), the ratio of the Roche-lobe 
radius and the orbital separation {Rl/ci) and the masses of both components (Mopt & Mx) of the three systems. The mean value and Icr Umits are 
shown with soUd and dotted Unes, respectively. 



of the compact object (see also Joss & Rappaport (1984) for a 
review on mass determinations in X-ray binaries). 

The ratio of the Roche-lobe radius and the orbital separation 
can be approximated by: 

— «i A-i-Blog^-i-Clog^^. (9) 
a 

The values of the constants A, B and C were determined by 
Rappaport & Joss (1983) to be: 

A ^ 0.398 - 0.026Q^ -i- 0.004i22 (10) 

B « -0.264 -I- 0.052Q^ - O.OISQ^ (11) 

C «i -0.023 - 0.005Q^ (12) 

where Q. is the ratio of the rotational frequency of the optical 
companion to the orbital frequency of the system; in case of syn- 
chronous rotation Q = 1. However, the timescale at which these 



systems are expected to synchronise is slightly longer than the 
timescale at which the orbit will circularise (for a detailed de- 
scription see e.g. Hut 1981). Therefore, these systems may still 
be in the process of synchronising the optical companion to the 
orbit, whereas their orbits have already become circular. The fact 
that the orbital periods of all three systems are decreasing, sug- 
gests that the donor stars are rotating slower than synchronous 
and that tidal forces are transferring orbital angular momentum 
to synchronise the system. 

It is possible to determine D. by measuring the projected ro- 
tation velocity Vrot sin i of the OB companion using the spectra at 
<p = 0.0. We use the grid of unified stellar atmosphere/wind mod- 
els of early-type supergiants computed by Lenorzer et al. (2004) 
using CMFGEN (Hillier & Miller 1998). First we select the lines 
that correspond to Unes included in the model. Subsequently, we 
select the model atmosphere that reproduces the observed line 
spectrum best. The normalised flux of the model is then scaled 
to yield exactly the observed EW. The models that correspond 
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Table 6. List of final values for the radial velocity amplitude (K„^i) and the results from the output of the Monte-Carlo simulations, i.e. the values 
of the inclination angle (i), the ratio of the rotational and orbital frequency (fl), the ratio of the Roche-lobe radius and the orbital separation (Ri^/a), 
the radius of the optical companion (i?opt), the orbital separation (a), and the masses of both components (Mppt & Mx) of the three systems. All 
errors are 1 cr. 



SMC X- 1 LMC X-4 Cen X-3 



^:opt(kms~') 20.2 ±1.1 35.1 ± 1.5 27.5 ± 2.3 

!(deg) 67+5 68+4 72+6 

a 0.91+[!ji 0.97+° !* 0.75 ±0.13 

RJa 0.61!;| 0.59!™ 0.63!««^ 

/fop.(Ro) 16.4!!? 7.8tl 12.1+0.5 

a(Ro) 26.6 ±0.8 13.3 ±0.3 19.r°l 

Mop, (Mo) 15.7:!-5 14.5+! i 20.2+!| 

Mx (Mq) i.06!°-{' l-25!°:i' 1-34!°:|^ 



best to our selection of lines are named "ARlIa", "ARIIII" and 
"09m-AR2m" for SMC X-1, LMC X-4 and Cen X-3, respec- 
tively by Lenorzer et al. (2004). Note that the names of these 
models correspond to a set of model parameters describing the 
model and not to the observational spectral type naming conven- 
tion. These models are subsequently convolved with a rotational 
broadening profile with a limb-darkening coefficient 0.6, as de- 
scribed by Gray (1992), to determine the value of Vrot sin ;. This 
results in 170 ± 30 km s'K 240 ± 25 km s"' and 200 ± 40 km s"^ 
for SMC X-1, LMC X-4 and Cen X-3, respectively. The se- 
lected lines and their corresponding best model are shown in 
Fig. 10. The rotational velocities required for a synchronous or- 
bit are 185 km s"', 250 km s"' and 255 km s"' for SMC X-1, 
LMC X-4 and Cen X-3, respectively. We will find below that 
these correspond to rotation rates consistent with, though per- 
haps slightly slower than, a synchronous orbit. 

Since for Roche-lobe overflow systems /3 > 0.9 (Avni & 
Bahcall 1975), we follow the approach of Rappaport & Joss 
(1983) and adopt that j6 is in the range 0.9 - 1.0. Thus, given 
a set of ^opt> P> «x sin i and 6^, we can determine by means of 
Monte-Carlo simulations a Icr confidence range for the values of 
Rh/a, i, Mopt andMx (see Rappaport & Joss 1983; Van Kerkwijk 
etal. 1995a). 

Since in these systems soft X-rays are absorbed by the ex- 
tended steUar wind of the optical companion, the eclipse lasts 
longer at low energies (up to ~ 5 keV), depending on the den- 
sity structure of the stellar wind of the optical companion (e.g. 
for 4U 1700-37 the eclipse at energies up to ~ 2 keV lasts al- 
most twice as long as at ~ 6 keV; see Haberl et al. 1994; Van der 
Meer et al. 2005). Therefore, we prefer 6^ determinations ob- 
tained from X-ray observations at high energies (see Table 1). 
For Cen X-3 these measurements are available. An accurate 
modelling of the X-ray light curve of multiple observations ob- 
tained with SAS-3/XTCA in the energy range 7.9 - 20 keV by 
Clark et al. (1988) results in a value of 6*6 = 32.9° ± 1.4°. They 
do not list an error on their value, so we use the standard de- 
viation in their 0^ distribution. For SMC X-1 and LMC X-4 
no detailed modelling has been performed and measurements 
are mainly reported for older X-ray missions. We use the range 
26° - 30.5° for SMC X-1 based on observations of Primini et al. 
(1976) (28.2° ± 0.9°, 2-6 keV, SAS-3), Bonnet-Bidaud & Van 
der KUs (1981) (29.9° ± 0.2°, 2-12 keV, COS-B) and Schreier 
etal (1972a) (29.1°±2.8°, 2-6 keV, fZ/wM). ForLMC X-4 we 
adopt the range 25°-29° based on observations of Li et al. (1978) 
(29.0° ± 2.5°, 6-12 keV, SAS-3), White (1978) (26.2° ± 1.1°, 
2-16 keV, ARIEL V) and Pietsch et al. (1985) (27.1° ± 1.0°, 
2-7 keV, EXOSAT). 



Using the defined input distributions for the Monte-Carlo 
simulations, we can determine the distributions for ;, Q., Ri^/a, 
Ropt, a, Mopt and Mx- All the results are listed in Table 6 and 
some of the corresponding distributions are shown in Fig 1 1 . 
The masses of the neutron stars become: Mx = 1.06^q [q Mq 
for SMC X-1, Mx = ^■'25^Vio for LMC X-4 and Mx = 
1.34+q}^ M© for Cen X-3 at a Icr confidence level. Compared 
to the mass determinations listed in Sect. 2 the error on the neu- 
tron star mass is reduced by at least a factor two. 

The masses and radii of the OB companion stars all lie in 
the range of 15-20 Mq and 8-16 Rq, respectively, in line with 
previous values (Van Kerkwijk et al. 1995a). Conti (1978) and 
Kaper (2001) show that for these stars a higher mass is ex- 
pected based on their spectral classification and conclude that 
these stars are undermassive for their luminosity. This may be 
due to the phase of mass transfer prior to the supernova forming 
the neutron star in the system. A detailed modelling of the optical 
spectra of the OB companions will be presented in a forthcoming 
paper. 

6. The neutron star mass distribution 

Abubekerov et al. (2004) present a new method to derive the 
radial velocity curves for HMXB systems in which the optical 
component is deformed due to the (partial) filling of its Roche- 
lobe. They show that these systems are aff"ected by gravitational 
darkening and by X-ray heating of the surface of the optical com- 
ponent. This can result in an underestimate of the radial velocity 
amplitude of the optical component and therefore in an under- 
estimate of the mass of the neutron star. Since this will mostly 
affect systems hosting a bright X-ray source, Roche-lobe over- 
flow systems will suffer most from this effect, i.e. the systems 
discussed in this paper. With our VLT/UVES observations it is 
possible for the first time to accurately determine the radial ve- 
locity amplitude of each absorption fine separately. We showed 
that for fines that vary in EW the radial velocity is indeed un- 
derestimated by several km s~^, consistent with the predictions 
of Abubekerov et al. (2004). By applying our carefuUy chosen 
selection criteria, we anticipate for these effects. 

We conclude that with our observations we have significantly 
improved the accuracy of the determination of the radial- velocity 
amplitude, and subsequently the determination of the neutron 
star mass in these three systems. Whereas some HMXB sys- 
tems have shown to host a neutron star with a mass significantly 
higher than 1.4 Mq, as is the case for Vela X-1 and possi- 
bly 4U 1700-37, the mass of SMC X-1 is low, 1.06+j; [J Mq. 
The masses of LMC X-4 and Cen X-3 are 1.251^ ji and 
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Fig. 12. Neutron Star (NS) masses in X-ray binaries, NS-NS binaries and NS-White Dwarf (WD) binaries obtained from Stairs (2004) and 
references therein. The indicated masses of SMC X-1, LMC X-4 and Cen X-3 are obtained from this study. The error bars correspond to 1 cr 
errors. This plot clearly suggests that neutron stars do not all have the same "canonical" mass. Note that although 4U 1700-37 most probably is a 
neutron star, it could be a black hole. 



1.34|J;pJ^ Mq, respectively. The mass of SMC X-1 is just above 
the minimum neutron star mass of ~ 1 Mq and significantly 
different from the mass of the neutron star in Vela X- 1 . We con- 
clude that the neutron stars in HMXBs have different masses, 
i.e. they do not all have the same "canonical" mass. We illus- 
trate our new mass derivations in Fig. 12, as part of the neutron 
star masses reported by Stairs (2004) and references therein for 
neutron stars in different types of systems. 

7. Discussion 

It remains to be explained why the mass of SMC X-1 is well 
below 1.28 Mq. The low mass may be the result of a different 
formation scenario, i.e. the electron-capture collapse of a degen- 
erate O -Ne -Mg core. Van den Heuvel (2004) argues that the 
generally low masses of neutron stars measured in binary radio 
pulsar systems may be due to a selection effect, as follows. Pfahl 
et al. (2002) noticed that there are two classes among the wide 
Be/X-ray binaries: (1) a substantial group with low orbital eccen- 
tricities, which indicates that their neutron stars received hardly 
any velocity kick in their formation events, and (2) a group with 
high orbital eccentricities, in which the neutron stars must have 
received a kick velocity of several hundreds of km s"' in their 
birth events. 

It was subsequently noticed (Van den Heuvel 2004) that the 
low orbital eccentricities of 5 out of the 7 known double neutron 
stars in the galactic disc indicate that the second-born neutron 
stars in these systems received hardly any kick velocity during 
their birth events and thus appear to belong to the same low- 
kick class of neutron stars as the ones in the low-eccentricity 
Be/X-ray binaries. These second-born neutron stars in the low- 



eccentricity double neutron star systems all appear to have low 
masses, in the range 1.18 - 1.36 Mq. This fits excellently with 
neutron-star formation by the electron-capture collapse of a de- 
generate O -Ne -Mg core, which is expected to form at the end of 
the evolution of stars that originated in the main-sequence mass- 
range 8 to about 13 ± 1 Mq (Miyaji et al. 1980; Podsiadlowski 
et al. 2005; Kitaura et al. 2005). Stars with larger masses de- 
velop at the end of their lives a collapsing iron core, surrounded 
by convective shells with O - and Si -burning. The violent con- 
vection in these shells may create large density inhomogeneities 
in the layers surrounding the proto-neutron-star formed by the 
collapsing iron core. This may lead to large anisotropics in 
the neutrino transport through these layers, which may cause 
the neutron star to be imparted with a "kick" velocity of some 
500 km s-HBurrows & Hayes 1996; Scheck et al. 2004). Indeed, 
young single radio pulsars have large space velocities (Gunn & 
Ostriker 1970) and their velocity distribution is very well repre- 
sented by a Maxwellian with a characteristic mean velocity of 
about 400 km s"' (Hobbs et al. 2005). 

In the light of these findings the low mass of the neutron star 
in SMC X-1 would be consistent with its formation by electron- 
capture collapse in a degenerate O-Ne-Mg core. This would 
imply a main-sequence progenitor mass < 14 Mq. Presently the 
companion of SMC X-1 has a mass of about 16 Mq. Allowing 
for some mass loss by stellar wind, its mass just after the mass 
transfer and the formation of the neutron star would have been 
about 18 Mq. 

With an explosive mass loss during the formation of the neu- 
tron star of about 1 Mq and a few solar masses stellar wind mass 
loss from the neutron-star progenitor, the initial system must 
have had a mass > 22 Mq (including the neutron-star mass). 
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Thus, a progenitor system of 13 Mo+9 Mq (or 14 Mq + 8 M©) 
would be consistent with the present system configuration. A po- 
tential problem with such a configuration is that conservation of 
mass and orbital angular momentum during mass transfer would 
lead to a fairly wide presupernova system, such that the present 
orbital period of ~ 3.9 days would be hard to understand, unless 
a large amount of orbital angular momentum has been lost with 
relatively little mass (at most a few solar masses) from the sys- 
tem. We thus conclude that the low mass of the neutron star in 
SMC X- 1 is consistent with formation by electron-capture col- 
lapse, provided that a relatively large amount of orbital angular 
momentum was lost from the system during the first phase of 
mass transfer. 

Acknowledgements. AvdM is supported by the Nederlandse Onderzoekschool 
voor Astronomie (NOVA). We would like to thank Rohied Mokiem for help- 
ing us with the grid of unified stellar atmosphere/wind models and Godelieve 
Hammerschlag-Hensberge for constructive discussions. The ESO Paranal staff 
is acknowledged for carrying out the sevice mode VLT/UVES observations. We 
are grateful to the anonymous referee for his/her constructive comments that 
helped to improve the quality of the paper. 

References 

Abubekerov, M. K., Antokhina, E. A., & Cherepashchuk, A. M. 2004, 

Astronomy Reports, 48, 89 
Ash, T. D. C, Reynolds, A. R, Roche, R, et al. 1999, MNRAS, 307, 357 
Aslanov, A. A. & Cherepashchuk, A. M. 1982, AZh, 59, 290 
Avni, Y. & Bahcall, J. N. 1975, ApJ, 197, 675 

Barziv, O., Kaper, L., Van Kerkwijk, M. H., Telling, J. H., & Van Paradijs, J. 

2001, A&A, 377, 925 
Bildsten, L., Chakrabarty, D., Chiu, J., et al. 1997, ApJS, 113, 367 
Blondm, J. M. 1994, ApJ, 435, 756 

Bonnet-Bidaud, J. M. & Van der Klis, M. 1981, A&A, 97, 134 
Brown, G. E. & Bethe, H. A. 1994, ApJ, 423, 659 
Burrows, A. 2000, Nature, 403, 727 

Burrows, A. & Hayes, J. 1996, Physical Review Letters, 76, 352 
Chevalier, C. & llovaisky, S. A. 1977, A&A, 59, L9 

Chodil, G., Mark, H., Rodrigues, R., et al. 1967, Physical Review Letters, 19, 
681 

Clark, G. W., Minato, J. R., & Mi, G. 1988, ApJ, 324, 974 
Clarkson, W. I., Charles, P A., Coe, M. J., et al. 2003, MNRAS, 339, 447 
Conti, P S. 1978, A&A, 63, 225 
Conti, P S. & Alschuler, W. R. 1971, ApJ, 170, 325 
Crampton, D., Hatchings, J. B., & Cowley, A. P 1985, ApJ, 299, 839 
Dekker, H., D'Odorico, S., Kaufer, A., Delabre, B., & Kotzlowski, H. 2000, in 
Proc. SPIE Vol. 4008, p. 534-545, Optical and IR Telescope Instrumentation 
and Detectors, Masanori lye; Alan F. Moorwood; Eds., 534-545 
FuUerton, A. W., Gies, D. R., & Bolton, C. T. 1996, ApJS, 103, 475 
Giacconi, R., Gursky, H., Kellogg, E., Schreier, E., & Tananbaum, H. 1971, ApJ, 
167, L67+ 

Giacconi, R., Murray, S., Gursky, H., et al. 1972, ApJ, 178, 281 
Gies, D. R., Bagnuolo, W. G., & Penny, L. R. 1997, ApJ, 479, 408 
Gray, D. F. 1992, The Observation and Analysis of Stellar Photospheres 
(The Observation and Analysis of Stellar Photospheres, by David F. Gray, 
pp. 470. ISBN 0521408687. Cambridge, UK: Cambridge University Press, 
June 1992.) 

Gunn, J. E. & Ostriker, J. P 1970, ApJ, 160, 979 

Haberl, F, Aoki, T, & Mavromatakis, F 1994, A&A, 288, 796 

Haensel, P, Zdunik, J. L., & Douchin, F 2002, A&A, 385, 301 

Heemskerk, M. H. M. & Van Paradijs, J. 1989, A&A, 223, 154 

Hilditch, R. W., Howarth, 1. D., & Harries, T. J. 2005, MNRAS, 357, 304 

Hillier, D. J. & Miller, D. L. 1998, ApJ, 496, 407 

Hiltner, W. A., Werner, J., & Osmer, P 1972, ApJ, 175, L19+ 

Hobbs, G., Lorimer, D. R., Lyne, A. G., & Kramer, M. 2005, MNRAS, 360, 974 

Hut,P 1981, A&A, 99, 126 

Hutchings, J. B., Cowley, A. P., Crampton, D., Van Paradus, J., & White, N. E. 

1979, ApJ, 229, 1079 
Hutchings, J. B., Crampton, D., & Cowley, A. P 1978, ApJ, 225, 548 
Jacoby, B. A., Hotan, A., Bailes, M., Ord, S., & Kulkami, S. R. 2005, ApJ, 629, 

L113 

Jones, C, Forman, W., Tananbaum, H., et al. 1973, ApJ, 181, L43-I- 
Joss, P C. & Rappaport, S. A. 1984, ARA&A, 22, 537 

Kaper, L. 2001, in ASSL Vol. 264: The Influence of Binaries on Stellar 
Population Studies, 125-+ 



Kaper, L., Hammerschlag-Hensberge, G., & Zuiderwijk, E. J. 1994, A&A, 289, 

846 

Kaper, L. & Van der Meer, A. 2005, ArXiv Astrophysics e-prints 
Kelley, R. L., Jemigan, J. G., Levine, A., Petro, L. D., & Rappaport, S. 1983, 
ApJ, 264, 568 

Kitaura, F. S., Janka, H. ., & Hillebrandt, W. 2005, ArXiv Astrophysics e-prints 
Krzeminski, W. 1974, ApJ, 192, L135 

Lang, F. L., Levine, A. M., Bautz, M., et al. 1981, ApL 246, L21 
Lattimer, J. M. & Prakash, M. 2004, Science, 304, 536 
Lennon, D. J., Dufton, P L., & Fitzsimmons, A. 1993, A&AS, 97, 559 
Lenorzer, A., Mokiem, M. R., de Koter, A., & Puis, J. 2004, A&A, 422, 275 
Levine, A., Rappaport, S., Deeter, J. E., Boynton, P. E., & Nagase, F. 1993, ApJ, 
410, 328 

Levine, A. M., Rappaport, S. A., & Zojcheski, G. 2000, ApJ, 541, 194 
Li, F, Rappaport, S., & Epstein, A. 1978, Nature, 271, 37 
Liller, W. 1973, ApJ, 184, L37+ 

Lutovinov, A., Revnivtsev, M., Gilfanov, M., et al. 2005, A&A, 444, 821 
Mason, K. O., Branduardi, G., & Sanford, P 1976, ApJ, 203, L29 
Miyaji, S., Nomoto, K., Yokoi, K., & Sugimoto, D. 1980, PASJ, 32, 303 
Mokiem, M. R., de Koter, A., Puis, J., et al. 2005, A&A, 441, 711 
Nagase, F, Corbet, R. H. D., Day, C. S. R., et al. 1992, ApJ, 396, 147 
Negueraela, 1., Smith, D. M., Reig, P., Chaty, S., & Torrejon, J. M. 2006, in 

ESA Special Publication, Vol. 604, ESA Special Publication, ed. A. Wilson, 

165-170 

Nice, D., Splaver, E., Stairs, 1., et al. 2005, ArXiv Astrophysics e-prints 
Paul, B., Raichur, H., & Mukherjee, U. 2005, ArXiv Astrophysics e-prints 
Petro, L. D. & Hiltaer, W. A. 1982, NASA STI/Recon Technical Report N, 85, 
19905 

Pfahl, E., Rappaport, S., Podsiadlowski, P, & Spruit, H. 2002, ApJ, 574, 364 
Pietsch, W., Voges, W., Pakull, M., & Staubert, R. 1985, Space Science Reviews, 
40, 371 

Podsiadlowski, R, Dewi, J. D. M., Lesaffre, R, et al. 2005, MNRAS, 361, 1243 
Priedhorsky, W. C. & Terrell, J. 1983, ApJ, 273, 709 
Primini, F, Clark, G. W., Lewin, W., et al. 1976, ApJ, 210, L71 
Quaintrell, H., Norton, A. J., Ash, T. D. C, et al. 2003, A&A, 401, 313 
Rappaport, S. A. & Joss, P. C. 1983, in Accretion-Driven Stellar X-ray Sources, 
1-39 

Reynolds, A. P, Bell, S. A., & Hilditch, R. W. 1992, MNRAS, 256, 631 
Reynolds, A. P, Hilditch, R. W., Bell, S. A., & Hill, G. 1993, MNRAS, 261, 337 
Reynolds, A. P., Owens, A., Kaper, L., Parmar, A. N., & Segreto, A. 1999, A&A, 

349, 873 
Sahade, J. 1959, PASR 71, 151 

Sahade, J. 1962, in Proceedings of a Symposium on Stellar Evolution, held in 
La Plata, November 7-11, 1960, La Plata: National University, Astronomical 
Observatory, 1962, edited by Sahade, J., p.l85, 185-+ 

Sanduleak, N. & Philip, A. G. D. 1976, lAU Circ, 3023, 1 

Savonije, G. J. 1978, A&A, 62, 317 

Savonije, J. 1983, in Accretion-Driven Stellar X-ray Sources, 343-366 
Scheck, L., Plewa, T, Janka, H.-T, Kifonidis, K., & Miiller, E. 2004, Physical 

Review Letters, 92, 01 1 103 
Schreier, E., Giacconi, R., Gursky, H., Kellogg, E., & Tananbaum, H. 1972a, 

ApJ, 178, L71 + 

Schreier, E., Levinson, R., Gursky, H., et al. 1972b, ApJ, 172, L79-H 
Schrijvers, C. & Telting, J. H. 1999, A&A, 342, 453 

Schrijvers, C, Telting, J. H., Aerts, C, Ruymaekers, E., & Henrichs, H. F. 1997, 

A&AS, 121, 343 

Srinivasan, G. 2001, in Black Holes in Binaries and Galactic Nuclei, 45—1- 

Stairs, I. H. 2004, Science, 304, 547 
Struve, O. 1937, ApJ, 85,41 

Thorsett, S. E. & Chakrabarty, D. 1999, ApJ, 512, 288 

Timmes, F X., Woosley, S. E., & Weaver, T. A. 1996, ApJ, 457, 834 

Tjemkes, S. A., Van Paradijs, J., & Zuiderwijk, E. J. 1986, A&A, 154, 77 

Val Baker, A. K. F, Norton, A. J., & Quaintrell, H. 2005, A&A, 441, 685 

Van den Heuvel, E. P J. 2004, in ESA SP-552: 5th INTEGRAL Workshop on 

the INTEGRAL Universe, 185-+ 
Van der Meer, A., Kaper, L., di Salvo, T, et al. 2005, A&A, 432, 999 
Van Kerkwijk, M. H., Van Paradijs, J., & Zuiderwijk, E. J. 1995a, A&A, 303, 

497 

Van Kerkwijk, M. H., Van Paradijs, J., Zuiderwijk, E. J., et al. 1995b, A&A, 303, 

483 

Van Paradijs, J. A., Hammerschlag-Hensberge, G., & Zuiderwijk, E. J. 1978, 

A&AS, 31, 189 

Vidal, N. v., Wickramsinghe, D. T, & Peterson, B. A. 1973, ApJ, 182, L77-I- 

White, N. E. 1978, Nature, 271, 38 

Wojdowski, R, Clark, G. W., Levine, A. M., Woo, J. W., & Zhang, S. N. 1998, 

ApJ, 502, 253 



A. van der Meer et al.: Determination of the mass of the neutron star in SMC X-1, LMC X-4 and Cen X-3 19 



List of Objects 

'SMC X-r on page 1 
'LMC X-4' on page 1 
'Cen X-3' on page 1 
'SMC X-r on page 1 
'LMC X-4' on page 1 
'Cen X-3' on page 1 
'SMC X-r on page 1 
'LMC X-4' on page 1 
'Cen X-3' on page 1 
'SMC X-r on page 1 
'LMC X-4' on page 1 
'Cen X-3' on page 1 
'SMC X-r on page 1 
'SMC X-r on page 1 
'LMC X-4' on page 1 
'Cen X-3' on page 1 
'VelaX-l' on page 1 
'J075 1 + 1807' on page 1 
'Vela X- 1 ' on page 1 
'SMC X-r on page 2 
'LMC X-4' on page 2 
'Cen X-3' on page 2 
'SMC X-r on page 2 
'LMC X-4' on page 2 
'Cen X-3' on page 2 
'SMC X-r on page 2 
'LMC X-4' on page 2 
'Cen X-3' on page 2 
'J1909-3744' on page 2 
'J075 1 + 1807' on page 2 
'Vela X-1' on page 2 
'SMC X-r on page 2 
'LMC X-4' on page 2 
'Cen X-3' on page 2 
'Vela X- 1 ' on page 2 
'4U 1538-52' on page 2 
'SMC X-r on page 2 
'LMC X-4' on page 2 
'Cen X-3' on page 2 
'4U 1700-37' on page 2 
'HD 153919' on page 2 
'Vela X-r on page 2 
'QV Nor' on page 2 
'4U 1538-52' on page 2 
'4U 1538-52' on page 2 
'SMC X-r on page 2 
'LMC X-4' on page 2 
'Cen X-3' on page 2 
'SMC X-r on page 3 
'LMC X-4' on page 3 
'Cen X-3' on page 3 
'SMC X-r on page 3 
'Sk 160' on page 3 
'SMC X-r on page 3 
'Sk 160' on page 3 
'SMC X-r on page 3 
'LMC X-4' on page 3 
'LMC X-4' on page 3 
'LMC X-4' on page 3 
'Cen X-3' on page 3 
'Cen X-3' on page 3 
'V779 Cen' on page 3 



'SMC X-r on page 4 
'LMC X-4' on page 4 
'Cen X-3' on page 4 
'SMC X-r on page 4 
'Cen X-3' on page 4 
'SMC X-r on page 5 
'LMC X-4' on page 5 
'Cen X-3' on page 5 
'SMC X-r on page 5 
'SMC X-r on page 6 
'LMC X-4' on page 6 
'Cen X-3' on page 6 
'SMC X-r on page 6 
'SMC X-r on page 7 
'V779 Cen' on page 7 
'SMC X-r on page 7 
'LMC X-4' on page 7 
'Cen X-3' on page 7 
'SMC X-r on page 7 
'LMC X-4' on page 7 
'Cen X-3' on page 7 
'SMC X-r on page 8 
'SMC X-r on page 8 
'SMC X-r on page 9 
'LMC X-4' on page 9 
'Cen X-3' on page 9 
'SMC X-r on page 10 
'SMC X-r on page 10 
'SMC X-r on page 10 
'LMC X-4' on page 10 
'Cen X-3' on page 10 
'SMC X-r on page 11 
'LMC X-4' on page 11 
'Cen X-3' on page 11 
'SMC X-r on page 11 
'SMC X-r on page 11 
'SMC X-r on page 11 
'V779 Cen' on page 11 
'Cen X-3' on page 11 
'SMC X-r on page 12 
'LMC X-4' on page 12 
'Cen X-3' on page 12 
'SMC X-r on page 12 
'LMC X-4' on page 12 
'Cen X-3' on page 12 
'SMC X-r on page 12 
'LMC X-4' on page 12 
'Cen X-3' on page 12 
'SMC X-r on page 12 
'LMC X-4' on page 12 
'Cen X-3' on page 12 
'SMC X-r on page 13 
'SMC X-r on page 13 
'LMC X-4' on page 13 
'Cen X-3' on page 13 
'SMC X-r on page 13 
'LMC X-4' on page 13 
'Cen X-3' on page 13 
'SMC X-r on page 13 
'LMC X-4' on page 13 
'Cen X-3' on page 13 
'Vela X-1' on page 13 
'SMC X-r on page 13 
'LMC X-4' on page 13 



20 



A. van der Meer et al.: Determination of the mass of the neutron star in SMC X-1, LMC X-4 and Cen X-3 



'Cen X-3' on page 13 
'SMC X-r on page 13 
'SMC X-r on page 14 
'LMC X-4' on page 14 
'Cen X-3' on page 14 
'SMC X-r on page 15 
'LMC X-4' on page 15 
'Cen X-3' on page 15 
'SMC X-r on page 16 
'LMC X-4' on page 16 
'Cen X-3' on page 16 
'SMC X-r on page 16 
'LMC X-4' on page 16 
'Cen X-3' on page 16 
'SMC X-r on page 16 
'LMC X-4' on page 16 
'Cen X-3' on page 16 
'SMC X-r on page 16 
'LMC X-4' on page 16 
'Cen X-3' on page 16 
'4U 1700-37' on page 16 
'Cen X-3' on page 16 
'SMC X-r on page 16 
'LMC X-4' on page 16 
'SMC X-r on page 16 
'LMC X-4' on page 16 
'SMC X-r on page 16 
'LMC X-4' on page 16 
'Cen X-3' on page 16 
'SMC X-r on page 16 
'LMC X-4' on page 16 
'Cen X-3' on page 16 
'SMC X-r on page 17 
'LMC X-4' on page 17 
'Cen X-3' on page 17 
'4U 1700-37' on page 17 
'SMC X-r on page 17 
'SMC X-r on page 17 
'SMC X-r on page 17 
'SMC X-r on page 17 
'SMC X-r on page 18 



